ACCRETION AND EMISSION PROCESSES IN AGN 



S. COLLIN 

Observatoire de Paris-Meudon, 92195 Meudon, France 
E-mail: suzy.collin@obspm.fr 

The UV-X continuum, the X-ray spectral features, and the variabihty in these 
bands provide powerful tools for studying the innermost regions of AGNs from 
which we gain an insight into the accretion process. In this chapter the discussion 
focusses on luminous AGN, i.e. Seyfert galaxies and quasars. The standard accre- 
tion disk model (a stationary geometrically thin disk) is described, and vertically 
averaged solutions for the radial structure are given. The emission of the standard 
disk is discussed using different approximations, and it is compared to the observa- 
tions. This leads to the conclusion that more complex models are required, such as 
the irradiated disk and the disk-corona models. The advantage of this last model is 
that it explains the overall UV-X spectral distribution. In the framework of these 
disk models, the profile, intensity, and variability properties of the X-ray iron line 
can be explained by reprocessing at the surface of the cold disk very close to the 
black hole (the "relativistic disk model"). An alternative possibility is discussed, 
where the UV-X continuum is produced by a quasi-spherical distribution of dense 
clouds surrounded by (or embedded in) a hot medium. In such a model the iron 
line profile could be due to Comptonization instead of relativistic effects. 



1 A Few Introductory Words 

In these lectures, I adopt the "black hole parajdigm" , in which AGN are fueled 
by accretion onto a massive black hole (BH)ll23. Besides the BH paradigm, 
there is another generally accepted model, i.e. that accretion takes place via 
a disk. A few initial comments should be made concerning this point. 

One tries to account for AGN luminosity by invoking accretion with an- 
gular momentum. If the flow is rotationally supported, it will form a disk. 
The presence of a disk is inferred from several observational facts. The most 
obvious one is the existence of a privileged direction: several active galaxies 
display cones of ionized gas, and collimated jets are commonly observed in 
radio-loud objects, as well as in a few radio-quiet ones. The "Unified Scheme"la 
implies the existence of a dusty torus. Double peaked line profiles observed in 
a small fraction of AGN are also believed to be the signature of disk emission. 
Finally, large gaseous disks have been observed in some cases, for instance 
in NGC 4258. However, one should keep in mind that all these evidences 
point towards the existence of a disk at large distances, i.e. at i? ^ 10'*i?G 
(i?G = GM/c^ is the gravitational radius, where M is the BH mass), while the 
bulk of the accretion luminosity is emitted by a very small region, at 10i?G- 
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Figure 1. Schematic view of an AGN, showing the dimensions of different regions in Rq 
(and in pc for M = 10* M0). The present chapter deals with the region located inside the 
"gravitationally unstable" part of the disk, closer than a few 10'^ -Rq from the BH. 



At such small distances the only evidence for an accretion disk is the profile 
of the Fe K line. Since this issue is important (as it implies the presence of a 
"cold disk" at only a few Rq in some cases), one should question this "proof, 
and examine whether there is not another possible interpretation. Moreover, 
there are some theoretical hints for quasi-spherical accretion. For these rea- 
sons I discuss quite extensively the "quasi-spherical" model, which docs not 
preclude the presence of angular momentum, although this is presently not a 
popular model. 

In order to show the physical scale of the region discussed in this chapter. 
Fig. 1^ displays a very schematic representation of the central flow in an AGN. 
We shall deal with the region located inside the "gravitationally unstable" 
part of the disk, closer than a few 10'^i?G (the distance where the broad lines 
are formed). There is presently no model proposed for the locally unstable 
region indicated in the figure, as the gravitational torques discussed in F. 
Combes' lectures concern the outermost parts of the nucleus, beyond 10 pc. 

I have chosen to discuss extensively the most simple models, and to men- 
tion only briefly the recent more sophisticated developments. My aim is also 
to stress the uncertainties underlying some widely accepted and fashionable 
models. Finally I should mention that I focus on radio-quiet objects, where 
accretion processes dominate the spectrum. 

The observational evidence is summarized in Sec. 2, and the basic physics 
of BHs are given in Sec. 3. In Sec. 4 a few results on the physical state of 
the medium emitting the bulk of the accretion power are deduced directly 
from the observations. In Sec. 5, the basis of accretion disk theory is re- 
called, emphasising some robust results concerning the emission spectrum, 
which are almost independent of the model. In Sec. 6, the standard (i.e. thin 
and steady) a-disks are discussed in their simplest (vertically averaged) form. 
Sec. 7 is devoted to the vertical structure. Sec. 8 tackles different models of 
"irradiated" disks, and their emission spectrum. Sec. 9 discusses the "cloud 
model", which is an alternative to these disks. In Sec. 10 the "thick" disk 
models envisioned for very low and very high accretion rates are briefly re- 
called. Sec. 11 reminds us that all this chapter deals with regions close to the 
BH, and that the accretion process is still not understood on larger scales. 



suzy-collin-mexico: submitted to World Scientific on February 1, 2008 2 



Figure 2. Typical AGN continuum of radio-loud and radio-quiet objects, and possible emis- 
sion mechanisms (some of them not mentioned in the text). Figure adapted from Koratkar 
& Blaes^'^, ©1999, Astr. Soc. of the Pacific, reproduced with permission of the editors. 



Figure 3. Mean line profile of the Fe K line for a sample of Sy 1 nuclei. The solid line is a 
double-Gaussian fit to the profile. From Nandra et al.*^. 



2 Summary of Observations 

I give here a very brief summary of some basic observational results which are 
crucial to constrain a model of the central engine|_For more details the reader 



2.1 The Broad-Band Spectrum 

Fig. ^ displays a typical broad-band continuunL from the radio to the 7- 
ray range, for radio-quiet and radio-loud AGNllij. The continuum in radio- 
quiet objects can be divided into three parts: an "IR-bump" , a "UV-bump" 
(creating a gap at about 1/im), and an X-ray power-law continuum. The 
"l/im gap" is clearly due to the superposition of two different components 
which dominate the IR and the UV emission, and fade at 1/im. Since l/im 
is the shortest wavelength at which hot dust can radiate, the IR-bump is 
attributed to hot dust radiation, at least in radio-quiet objects. 

The continuum produced directly by accretion processes is restricted to 
the UV-X range (though the mid-IR continuum could be partly reprocessed 
nuclear radiation). From 1 to 10 eV the spectral index a, defined as oc 
j/^", lies in the range from —0.3 (for luminous quasars) to +0.5 (for Seyfert 
galaxies and low luminosity quasars). The soft X-ray continuum below 1 keV 
seems to be simply the extension of the UV with a spectral index close to 
1.5. The continuum from 1 eV to about 100 cV is known as the Big Blue 
Bump (BBB), with a peak at 10 eV (see Fig. 5 in H. Netzer's chapter). A 
noticeable fact concerning the UV continuum is the almost complete absence 
of a Lyman discontinuity (a weak absorption edge is sometimes observed but it 
is most probably due to the Broad Line Region or to intervening intergalactic 
clouds with a redshift close to the emission one). I will not discuss here the 
"Warm Absorber" , which imprints features in the soft X-ray continuum, as it 
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is covered in H. Netzer's lecture notes. 

Above 1 keV the average spectral index is about 0.7, so that in the ma- 
jority of objects the continuum below 1 keV shows a "soft X-ray excess" when 
compared to the extrapolation of the continuum in the 1-10 keV range. Recent 
observations performed with Chandra show that the strength of the excess has 
probably been overestimated in the past. Above 1 keV the spectrum can be 
decomposed into an underlying power-law with a spectral index (in photon 
number) T ~ 1.9 and, superimposed on it, a broad emission line peaking at 
6.4 keV, identified with the Ka line of weakly ionized iron_and a "hump" 
peakirig_at about 30 keV (for a review see Mushotzky et alEj). It has been 
shownll23 that the strength of the hump, measured by a parameter R, corre- 
lates with r. The power-law continuum has a cut-off with an e-folding energy 
of 50 to 300 keV. 

The Fe K line displaita an asymmetrical profile with an extended red wing 
and a smaller blue wingEJ (see Fig. ||). The presence_Qf the broad red wing 
in the average spectrum has been recently questionedE£j. However in the two 
best studied objects, MCG-6-30-15 and NGC 3516, the red wing extends down 



13J 



to 4 keV. The strength of the line seems also to correlate with T 

Finally, another important result is that the BBB luminosity is larger 
than the hard X-ray luminosity, in particular in luminous quasars. 



2.2 Variability Properties 

Since B. Peterson focussed mainly on the optical-UV properties, it is necessary 
to recall briefly some aspects of X-ray variability relevant to this chapter. 

A well established result is that the characteristic time variability in the 
UV scales with luminosity: in Seyfert galaxies it is typically of the order of 
a few days, while it is several months in high luminosity quasars. The X-ray 
flux varies more rapidly than the UV, with typical time-scales of hours for 
Seyfert galaxies. 

The conclusions one can draw from correlations between continuum bands 
are complex and sometimes controversial. For instance, NGC 7469 shows 
variations in the UV and X-rays with similar large amplitudes but with a 
time-delay of the order of 2 days, the UV leading the X-rays at maxima. 
The_X-ray fiux displays in addition short term variations not seen in the 
UVl3. Moreover, the X-ray spectral index seems to be correlated with the 
UV In NGC 3516 the delay between UV and X-bands is larger than 2 

dayscJ. In_NGC 4051 no optical variations were observed during strong X-ray 
variationsEEI. The UV variations drive the optical variatimis with a time-delay 
< 0.2 days in NGC 41510, < 0.15 days in NGC 35160, and - 1 day (but 
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increasing with wavelengthEii'lli) in NGC 7469. Soft X-ray variatiopa are 
generally larger than hard X-ray ones, but this is not always the caseE^I. _La 
NGC 3516, they were strongly correlated with no measurable lag (< 0.15d)c2l. 

The centroid energy, the intensity, and the equivalent width of the Fe Ka 
line, are variable on short time-scales (down to ks). Correlations between 
the continuum and Fe K have been intensively looked for in a few objects 
with ASCA, RXTE, and_.now with Chandra, revealing a complex behavior 
(for a review see Nandrao). ASCA has observed rapid (1 to 10 ks) variations 
of Fe K iijL_|Qmfile and intensity in MCG-6-30-15O. This result has been 
questionedEj0 by recent observations of MCG-6-30-15 and NGC 5548 with 
RXTE, which show that the line flux is constant over time-scales of 50 to 
500 ks, while the profile is variable and the underlying continuum displays 
large flux and spectrai-szariations! The same result seems to be observed in 



NGC 3516. ReynoldalHi excludes a time-delay of 0.5 to 50 ks between the 
line and the continuum, and suggests that the line flux remains constant on 
time-scales between 0.5 and 500 ks in MCG-6-30-15. An analysis of a sample 



of Seyfert nuclei from the ASCA archival^ finds that in most cases changes in 



the line do not appear to track changes in the continuum. Unfortunately, the 
new results found with Chandra are not yet published at the time of writing 
these lectures. 

We summarize here the most important points that any model should be 
able to explain: 

• Spectral features: 

— the overall shape of the continuum: the BBB (but does it include or 
not a soft X-ray excess?), the absence of the Lyman discontinuity, 
the UV/X ratio, the X-ray hump, the power-law, the break at ~ 100 
keV, the Fe K line; 

— the correlation between R and F; 

— the very broad red wing of the iron line (at least in a few objects), 
and its centroid energy. 

• Variability properties: 

— the rapid variations (< ks) of the 2-lOkeV continuum; 

— the correlation between (and only between) long-term variations of 
UV and X-ray fluxes; 

— the very short (or absent) time-lag between optical and UV varia- 
tions. 
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And perhaps a few other, but not so well established, variability properties: 
the slow variation of R; the rapid changes of the Fe K profile; and the absence 
of correlation of the Fe K flux with the underlying continuum. 

Several other issues are not mentioned here, including the polarization 
properties, in particular, the abrupt increase of polarization level observed at 
750 A in quasars. 



3 Black Holes, Accretion Rates, and Luminosities 

It is beyond the scope oLtlupfhapter to discuss BH physics, for which excellent 



text books already exist 121 'El. For the purpose of this chapter it is only neces- 
sary to know that there are two types of BHs: non rotating or Schwarzschild 
BHs; and rotating or Kerr BHs. The innermost stable orbit of a Schwarzschild 
BH is equal to 6Rg, while it depends on the angular momentum, J, for a Kerr 
BH. Due to photon capture the maximum value of cJ/GM^ = 0.998 and not 
unity, and it corresponds to a radius of the last stable orbit equal to 1.24i?G- 
Inside the innermost stable orbit, a particle falls almost radially onto the BH. 
The binding energy of the last stable orbit determines the efhciency of energy 
conversion. For Schwarzschild and Kerr BHs this is ?/ — 0.057 and rj ~ 0.30 
respectively. 

Another important parameter is the "gravitational redshift" . This is due 
to the increase of proper-time with decreasing distance from the BH. Since 
r = R/Rq can be smaller for Kerr BHs, photons orbiting around a Kerr BH 
can be more redshifted than around a Schwarzschild BH. 

The Eddington luminosity, -f/Eddj describes the maximum available power 
due to accretion. It corresponds to the equality between the gravitational force 
exerted on protons and the radiative force. In a fully ionized gas, the opacity 
is due to Thomson scattering, and this leads to the well known expression 

^Edd — , (J-j 

where ctt is the Thomson cross section and /ie is the unit mass per electron. 
In units appropriate for quasars, this can be written as 

LEdd = 1.5 X lO^^Mg erg s'^ , (2) 

where Afg is the mass expressed in 10* Mq. Note that this expression is valid 
in principle only for spherical accretion of a fully ionized medium in steady 
state. 

The corresponding accretion rate is 

-^Edd 



MEdd = - 2.8 Mi 

c rj 
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Figure 4. T vs. S for an optically thin medium exposed to a mixture of power-law contin- 
uum (of index 0.8 in Fi/) and a 10 eV BB, for increasing fractions of the BB intensity from 
top to bottom. Prom Nayakshin & Kallman^®. 



We will often use the "Eddington ratio" m = M/M-^dd = L/L^dd- Some 
people define instead m = M / L^dd, which differs by rj from our expression. 

B. Peterson showed in his lectures that BHs in Seyfert nuclei have a mass 
range of 10^ to 10^ M©, with Eddington ratios of about 0.1. 

4 What Does the BBB Tell us About the Emitting Medium? 

The shape of the BBB implies thermal emission from an optically thick 
medium radiating as a black body (BB) at temperatures from a few 10'*K 
(optical emission) to a few lO^K (EUV emission). The medium can however 
be optically thin at frequencies larger than the peak of the local Planck curve. 

If the medium radiates like a BB, the size of the emitting region, i?BBBi 
is given by 

Lbbb = Min(l,T) a 47ri?|BB (^T^s > (4) 

where Lbbb is the luminosity of the BBB, r is the optical thickness, a is a 
factor of the order of unity that depends on the geometry of the emitting 
region (in the case of a disk a = 0.5), a is the Stefan constant, and Teg is the 
effective temperature. Since r is larger than unity, using the definition of Rq 
one gets 

-RbBB „„ 1/2 ^-2 ,T-l/2 /-^BBB 

rBBB = 30 a / Tg Mg , (5) 

where T5 is the effective temperature in units of lO^K. 

One can deduce that the density of the emitting medium is high. Let 
us assume that the UV emitting medium "sees" the X-ray source, and that 
the surface layers are heated by this source. The temperature T of an X- 
ray heated gas, assuming a spherical geometry, is governed by a parameter S 
defined as 



1 Lx 

eCkT A-kB? ' ^ ' 
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where Lx is the X-ray luminosity. Note that S is closely related to the ion- 
ization parameter U used by H. Netzer in his lectures. The S vs. |-T curve 
depends on the incident spectrum, as first shown by Krolik et alE^I. This 
curve is characterized by two thermally stable solutions: a cool phase up to 
about 2 X 10^ K, where the energy balance is due to atomic processes; and 
a hot phase close to the "Compton temperature" Tcomp, that corresponds to 
the balance between Compton heating and Compton cooling. As an illustra- 
tion, Fig. ^ displays the S vs. T curves for an optically thin gas illuminated 
by a mixture of power-law continuum (index 0.8 for F^) and a 10 eV BB, for 
increasing fractions of the BB intensity. One can derive that the fraction of 
soft photons decreases with an increasing Tcomp- For instance, the composite 
Laor et alE3 AGN continuum has Tcomp = 2.5 x 10^ K, corresponding to the 
lower curve on the figure. We shall come back to the S vs. T curve in the 
following sections. For the moment, we only need to know that whatever the 
fraction of BB emission, a temperature of lO^K is obtained for S ~ 1 — 10. 
Since S must be of the order of a few units, using Eq. (Q) one gets: 

ne ~ 3 X 10^6 Afg-i ( ( ' cm-3 . (7) 



^Edd 



\10Rg J 



Of course, things are more complicated if the medium is stratified and heated 
by non radiative processes, but this result remains basically correct. 

The suggestion that this thick, dense, and relatively cold medium emitting 
the BBB could—he identified with an accretion disk was majie—for the first 



time by Shieldslli3 for 3C 273, followed by subsequent papersO'lIj. There are 
however other alternative explanations which will be discussed at the end of 
the chapter. We now concentrate on the physics of "standard accretion disks" 
and their emission properties. 



5 Thin Steady Accretion Disks: the Simplest Approximation 

First we can show that if the disk radiates like a BB at ~ 10^ K, it is ge- 
ometrically thin, i.e. its scale height H is much smaller than its radius R. 
The vertical component of the gravitational acceleration at an altitude z is 
GMz/R^, so H ^ Rcs/Vc/,, where Cs is the sound velocity and the orbital 
velocity. Although the temperature inside the disk can be (and actually is) 
larger than lO^K, at 10-100i?G Cs ^ V^, which equals a fraction of the light 
velocity, so H/R ^ 1. 

Since the disk is geometrically thin, the problem is reduced to two di- 
mensions. The gas velocity has only two components, and Vr, the radial 
component. If the BH mass dominates the mass of the disk and the central 
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stellar cluster (which is the case for R < \Q'^Rq)^ is Keplerian, so the 
orbital frequency f^K = {GM/R^y/'^, and one deduces the basic equation 

H^Cs^K , (8) 

which is a simplified form of the hydrostatic equilibrium equation. 

At this stage we can determine the dynamical time tdyn, i-e. the time it 
takes for a perturbation to cross the disk vertically. According to Eq. (^), 

_ ^ _ 1 

idyn — — — t;— , (9) 
Cs "K 

which written in appropriate units gives 

tdyn 500M8r3/2 s. (10) 

This is the shortest time on which one can expect a change of structure to 
occur at a given radius. 

The surface density, E, and Vr are governed by the continuity equation 
and the equation of motion: 

r)y 1 8 

where Q is the torque between two adjacent rings. 

To simplify the discussion, we consider a stationary disk. This is the case 
if the disk is fueled at a constant rate at its outer radius during the time it 
takes for the gas to reach the BH (the "viscous time"), and if there is no 
matter falling radially onto the disk at smaller radii. Instabilities can however 
suppress the stationarity. Under those conditions Eq. ( pT| ) becomes 

M = 2nRT,Vr , (13) 

where the accretion rate M is constant along the radius; and Eq. ( p^ ) becomes 

g = R'^VIkM + constant , (14) 

which describes the constancy of the angular momentum flux. From Eq. ( p^ 
it is interesting to note that the thin disk assumption (i.e. Vr <C V^) implies 
a large value of the surface density, 

^^lO^^l^^ gcm-2. (15) 
ijVr r 

Since matter orbits in circles it requires a dissipation mechanism to lose 
energy and move towards the BH. The question is how to ensure accretion, or 
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Figure 5. Disk spectrum, assuming that the disk radiates locally like a BB. Boundary 
corrections were not taken into account. The index 2 and the exponential cut-ofT correspond 
to the Rayleigh-Jeans and to the Wien part of the spectrum respectively. 



in other words: what is the torque G? Microscopic viscosity is far too small, 
so one has to seek for a more efficient mechanism. Let us assume for the 
moment that there is a local unspecified viscosity. 

The kinematic viscosity v is related to the torque through the viscous 
stress tensor, giving the rate at which angular momentum is transported. For 
a thin disk, the only interesting coordinate of the stress tensor is 

where p is the density. The torque is thus equal to 

g = j Rd(i) J RTr4,dz = 2ttR^vY.^^ . (17) 

Note that since dVlYi/dR is negative, the torque is negative, and angular 
momentum is transported outwards, as required for accretion to take place. 

Finally, one can compute the rate of energy dissipated in the disk per unit 
surface, D{R). According to our assumption of a local viscosity, the potential 
energy is transformed into heat, so we would get 2ttD{R) — GMM / R^ . How- 
ever, we also have to take into account the conservation of angular momentum 
and the fact that the torque cancels at the radius of the last stable orbit. 
This leadSpto an extra factor (3/2)(l— y^Rin/R) (for a detailed demonstration 
see KrolikEa). Finally one gets 



where f{R) = (1 — y/Rin/R)- In a fully relativistic treatment the expression 
of /(i?) would be more complicated. It is important to realize that the ap- 
proximate R^^ dependence of the dissipation leads to a luminosity roughly 
proportional to 1/R. Therefore, the bulk of the gravitational energy is re- 
leased close to the BH, at a few Rq, but this does not mean that the bulk of 
the radiation necessarily comes from this region (cf . the "cloud model" ) . 

If the disk radiates like a BB, the temperature Tbb is equal to the effective 
temperature T^s- Since the energy is radiated by both faces of the disk, 
Tcs — (^/2(7)^/'*, which translates into 

Teft - 6 X 10''f]-^/'^f{R)^/^m^/^M-^^\-^/^ K, (19) 



suzy-collin-mexico: submitted to World Scientific on February 1, 2008 10 



where m is defined by Eq. (|^). Integrating the Planck function over the 
radius and neglecting the factor /(i?), one obtains the spectrum shown in 
Fig. 1^, where Tin and Tout are the temperatures at the inner and outer radii. 
According to the Planck law, at a frequency v — kT^/h, the regions with 
T < Ty, i.e. R > R{T,y), contribute negligibly to the emission. On the 
other hand, in the inner regions where T > T^, the spectrum at the frequency 
1/ is given by the Rayleigh- Jeans expression = 2v^c^^kTcs. Thus, T^ cx 
/flou?^ R^'^dR, and finally the emission at a frequency v is dominated by 
material at the radius R{T^). 

This spectrum should be corrected for relativistic effects close to the BH, 
in particular in the case of a Kerr BH, where the inner radius is small. The 
effect is to make the spectrum dependent on the inclination angle, as L^, is 
relativistically boosted towards high frequencies for edge on disks (cf. Fig. [T^ ). 

The fact that the bulk of the emission at a given wavelength is produced 
at the radius Ri, has important consequences for the variations of the different 
spectral bands. Since Bi,(T) peaks at A = 0.51/T(K) cm, one deduces from 



Eq. (19), that the emission at A^, expressed in microns, is produced mainly 
at a distance 

r ~ 650 XyS'^'^f{Rf'^m^'^M^^'^ , (20) 

or 

R ~ l.2\y^j^-^/^f{RY'^m^'^Ml'^ It-d. (21) 

From this we infer that the optical to near-UV emission corresponds to regions 
located at large distances from the BH. Comparing this result with the small 
time-delays between the UV and optical light curves (less than a fraction of a 
day), we see that it is marginally compatible with a causal link propagating 
at the speed of light. It is however incompatible with the dynamical time, 
according to Eqs. ([l^) and (po|). But before we accept this conclusion we 
should question whether it is legitimate to assume that the disk is radiating 
locally like a BB, and therefore we have to obtain information on the optical 
thickness of the disk, and more generally on its physical properties. In this 
aim we have to specify the viscosity mechanism. 

6 A Further Step: the a-Prescription, Vertically Averaged 
Models 

6.1 Disk Structure 

It is widely accepted that the mechanism for transportation of angular mo- 
mentum can be identified with "turbulent viscosity" according to the "a- 
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prescription" proposed bK Shakura & Sunyaevti3 (for a detailed discussion of 
a-disks see Frank et al.H). In this prescription it is assumed that the sizes 
of the turbulent eddies are smaller than the thickness of the disk and that the 
turbulence is subsonic, hence 

V = aCsH , (22) 

where a ^ 1 . 

The a-prescription has proven to give satisfying results when applied to 
disks in cataclysmic variables, where a relatively small value (< 0.1) of a 
is generally reauired. A magnetohydrodynamki^ instability pointed out by 
Chandrasekharll3 is able to provide this viscositytEl. Other prescriptions_hased 
on purely hydrodynamic (shear) instabilities are also under debatcE££l, in 

SXk- The 
foUow- 



which the corresponding value of a can be larger and vary^ 
viscosity can be parameterized with a fixed Reynolds numbc 



ing an approach suggested already by Lynden-Bell & PringicET Finally the 
viscosity might also be a non-local mechanism, for instance in the case of an 
organized magnetic field anchored in the disk. 
From Eqs. ( p^ ) and (p^, one finds 

K = ^«c,| , (23) 

which shows that the radial velocity is strongly subsonic, and a fortiori 
strongly sub-keplerian, as expected. 

The a-prescription is equivalent to the previous expression if and only if 
the disk is not self-gravitating. Eqs. (|l^) and (^2|) give 

Tr^ = ^apCsHQK . (24) 

Using Cs — HilKi, one gets 

Tr^^^aP , (25) 

where P is the sum of the radiative and gas pressure. 

If the ratio of the vertical component of the disk gravity to that of the 
central gravity, C = 47rG/5/n^ is of the order of unity, thenE3: 

Cs ^ HflKil + 0'^^ and Tr^ = ^aP{l + Cy^^^ ■ (26) 

Actually there is a debate about the value of P one should take in the stress 
tensor: is it the total pressure or only the gas pressure? For the moment we 
adopt the use of the total pressure. 
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Because it is a matter of confusion, it is necessary to stress here that 
T ~ Tmid-pianc and not Tcff. If the diffusion approximation holds between 
the two temperatures, then T ~ r^/'*Tofr, which imphes that T ^ Tcff in the 
optically thick region. 

All these expressions are vertically averaged. Since the disk is geometri- 
cally thin, these averaged quantities can be used to decouple the radial and 
vertical structures. Additional equations describe: 



The gas pressure: 

gas — 

where is the mean mass per particle. 



(27) 



The radiative pressure, which in the 2-stream approximation is given by 



,y ^rad 

^rad ^ 



1 1 



(28) 



where tr is the Rosseland optical thickness of the disk. 

• The flux radiated by each face of the disk@ 

^rad = o 1 I 1 , 1 ' (29) 

where rp is the Planck optical depth. Note that, in the optically thick 
case, Prad tcuds to its LTE value and Piad ~ (8/3)crT"'/rR, which is the 
well-known diffusion approximation. In the optically thin case, _Fiad ~ 
4crrVp. 

• The energy balance between the local viscous dissipation and the radia- 
tive cooling, 

F,ad = i?(i?)/2 . (30) 

These equations, when solved self-consistently with the appropriate func- 
tions for kr,(p, T), Kp(p, T) and np{p, T), yield the radial disk structure in the 
non self-gravitating region of the disk for given values of M, M , and a. One 
finds that the disk is divided into several regions, according to the dominating 
processes: radiation or gas pressure; Thomson, Kramers, or more complicated 
opacity; and self-gravity or no self-gravity. 

To show the importance of the opacity. Fig. ^ displays the radial depen- 
dence of several physical quantities as function of the radius for an a-disk 
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Figure 6. From left to right and top to bottom, density, surface density, Rosseland mean 
optical thickness, scale height, equatorial temperature and radiation to gas pressure ratio 
as functions of the radius for a vertically averaged a-disk. The thick solid line corresponds 
to a 2D model with realistic opacities. Courtesy of J.-M. Hure. 



Figure 7. Transition from radiation to gas pressure regime for a vertically averaged o-disk, 
for M = IO^Mq. The figure displays m versus r for different values of a. Courtesy of J.-M. 
Hure. 



Figure 8. Transition to the self-gravity regime for a vertically averaged o-disk. The figure 
displays Rgg/RQ versus m for different values of M. All models except one are computed 
for a = 0.1. The curves are labelled with logAf in units of Mq. 



in the case of a typical AGN. In these computations self-gravity has been 
neglected in order to focus on the other parameters. Different approxima- 
tions have been compared to a 2D computation with realistic opacities. We 
see that in the inner region {R < 10^^ cm, i.e. r < lO'^iJc), the opacity is 
better approximated by Thomson scattering, while in the outer region it is 
better approximated by Kramers' opacity. Note that the radiation pressure 
dominates in the Thomson regime (/3 = Pgas/frad)- 

Fig. shows that the transition between the radiation and gas pressure 
dominated regimes occurs at a larger radius for a larger accretion rate in 
Eddington units. This radius also increases with a decreasing viscosity pa- 
rameter a, which corresponds to an increasing density. For a typical AGN 
(m = 0.1 and M = IO^Mq) and a = 0.1, R ^ 2000i?G- Note that it occurs 
in the same region as the transition to the self-gravity regime (this is due to 
the rapid increase of density when the disk gravity begins to dominate). 

Finally the disk becomes self-gravitating at a radius Rgg corresponding 
to C ^ 1. As shown in Fig. ^, the transition occurs at i? ~ 100 — 1000i?G- 
It is located further away for a larger value of a, owing to the corresponding 
smaller value of the density. 

In summary the non self-gravitating part of the disk is generally reduced 
to the Thomson opacity-radiation pressure dominated regime. Neglecting the 
boundary function f{R) the different parameters vary as: 
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• H = constant 



• S and T (X R^/^ , 

• T cx i?-3/8 ^ 

• Pgas/Prad « i?21/8 . 

These relations explain why the transition from the gas pressure to the radi- 
ation pressure regime, and the transition to self-gravity, are so rapid (see the 
strong variation of ( and Pgas/-Prad as functions of the radius on Fig. ^. It is 
also interesting to note that the gas temperature decreases less rapidly than 
Teff (which is proportional to R~^/^). The scale height to radius ratio can be 
expressed in this region as 

which shows that the "geometrically thin" assumption does not hold for large 
to; and the Thomson thickness is given by 



A few important time-scales relevant to the external regions which emit 
at optical wavelengths, where self-gravity is not negligible, are defined below. 
The "viscous time" ivisc, basically the time it takes for matter to move to the 
center, can be expressed as (using in particular Eq. (p6|)): 

^^'^^ ^ Vr ^ a[R) fWTT? ■ ^^^^ 

The "thermal time" tth, i-e. the time it takes to radiate the thermal energy, 
can be expressed as 



PH yi + c 

The expression for the dynamical time tdyn has already been given in Eq. (^ . 
Correcting for self-gravity, it becomes 
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Figure 9. Viscously (left) and thermally (right) unstable regions for a vertically averaged 
«-disk, with a = 1, M = 3 X IO'^Mq, and m=0.03. The figure displays T (in K) and H (in 
cm) versus r. The solid line represent the unstable regions. From Hure*^. 



It is easy to see from these expressions that 

^visc '3> tth > tdyn ■ (36) 

For instance in the inner Thomson opacity - radiation pressure dominated 
region ^visc is given by 

= 10 a-'^Msm-^r^/^ s. (37) 

6.1.1 Stability of the Disk 

There are 3 main instabihties acting in the disk: 
1. the thermal instabihty, which occurs when 



d\nT J „ \ d\nT y 
2. the viscous instabihty, which occurs when 

/ SIuM' 



ainS 



< ; (39) 

R 



3. the gravitational instability, which occurs whenLHZ'^ 



Q-%<i. (40) 

In the self-gravitating regime, whereS p cx and S oc i?^, the disk 
shrinks, rapidly becomes strongly gravitationally unstable, and unstable frag- 
ments with a typical volume begin to collapse. Consequently, for R > Rgg 
the a viscosity does not hold. We stress that this gravitational instability 
is local, and does not require the disk mass to be of the order of the BH 
mass. The instability corresponds simply to the tidal limit 47rl] > MH/R^, 
or equivalently to Afdisk ^ MH/R, hence Afdisk ^ M. 

It can be shown that the thermal and the viscous instabilities are strongly 
linked, so they generally occur in the same region of the disk (see Fig. We 
see that the disk is both viscously and thermally unstable in a large fraction 
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of the inner region. This combined instabiUty should operate on a thermal 
time-scale, which is quite small in the inner regions (cf. Eqs. (|o|) and (jsj)). 
It should therefore induce rapid variations of EUV emission. Moreover as 
a-disks are also gravitationally unstable in the outer region, one deduces that 
in AGN, a-disks are unstable almost everywhere! 

6.1.2 "13-Disks" 

In the case where the stress tensor is^aaumed to be proportional to the gas 
pressure and not to the total pressureli3'La, the disks are known as "/3-disks" . 
Strictly speaking they should be called "a — /3-disks" ! They are denser than a- 
disks, and dominated by gas pressure even in the inner regions. Consequently 
they are not subject to the viscous/thermal instability, but to another instabil- 
ity due to ionization of hydrogen, which occurs in regions where the mid-plane 
temperature is of the order of lO^K. This instability induces a "limit-cycle" 
behavior, which is thought to explain the periodic variations of dwarf novae. 
In AGN it should operate on a much larger time-scale (10^ yj4Ti and could 
correspond to variations of 4 orders of magnitude in luminosity 



120 



6.2 Where a Fundamental Ingredient, Compton Scattering, Appears 

The solution of the vertically averaged disk showed that Thomson opacity 
dominates over free-free opacity in the inner regions. Diffusions are therefore 
important, and modify the equilibrium spectrum. Radiation is thermalized 
for a relatively small absorption optical thickness, when the effective optical 
thickness Teg = ^/T^hJJah^+^^ V^absTx > 1, where Tabs and tt are the 
absorption and diffusion optical thicknesses respectively. This_accurs in the 



inner regions. In this case the spectral distribution is given bytil 



/. = 2i3.J^«i?, . (41) 
For instance, if absorption is dominated by free-free processes, 

Kabs = 3.69 X 10^ N,N,Z^gs^^ (^1 - exp ' (42) 



and, using Eqs. (|4l|) and (|42[), one deduces that the intensity is given by 

cx i^T^/^ (43) 
in the Rayleigh- Jeans regime, where hv/kT < 1; or 
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in the Wien regime, where hv/kT > 1. The spectrum emitted by the inner- 
most regions wiU therefore differ from a local BB. 

However it is not sufficient to take into account Thomson scattering, as 
one can easily show that Compton scattering is also important. Since the gas 
temperature is high and we are interested in soft X-ray photons, the collisions 
between electrons and photons are not elastic, and they lose or gain energy 
with each scattering. The "Compton parameter" y, which is the total relative 
change of energy when a photon travels through the medium, is equal to 
the mean number of diffusions times the relative change of energy at each 
scattering. The mean number of diffusions, Ns = Max(TT,T,f) = t^- In the 
case of a thermal non relativistic distrihulion of electrons, the relative change 
of energy at each scattering is given bjiiiil 



and thus one gets 



A/ii^ hv AkT , , 

-i— = 2+ 2 ' 45 



AkT — hv ^ , 

v = —r^^T ■ (46) 



In the inner regions of the disk, tt ^ lO"* and T ~ lO^K, so y can easily 
be larger than unity. If free-free dominates, then T(y) and hence y increases 
with decreasing frequency, and there is a frequency Vq below which y > \. 
If hvo/kT <^ 1, inverse Compton processes will drive soft photons to higher 
frequencies, and there will be an accumulation. joLphotons aX v ^ AkT/h, at 



the expense of the low frequency populationO^liHiJ. Another consequence of 
Compton diffusion is that the profile of a line will be asymmetrically broad- 
ened, as it is dominated by Inverse Compton diffusions for vunc < AkT /h, 
and by direct Compton diffusions for j^iino > AkT/h (i.e. down-scatterings 
leading to a broad red wing) . The computation of the spectrum in this case is 
performed with transfer codes using the "Kompaneets diffusion equation" , or 
with Monte Carlo codes. The effect of Compton scatterings on the spectral 
distribution of the disk is illustrated in Fig. |ll|, and at the end of the chapter 
we will discuss the effect on the profile of the iron line. 

Another problem which has to be taken into arr-mint is the existence of an 
atmosphere above the disk. As a first approachE3'lHj, the disk was replaced 
by low gravity stellar atmospheres of different Toff. However, since gravity is 
constant in a stellar atmosphere, while it increases with height in disks, this 
method is not valid, and led to the erroneous result that the spectrum should 
always show a strong Lyman discontinuity in absorption. 
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Figure 10. Temperature and density structure of an a-dislj for a typical AGN (M = IO^Mq, 
M = O.IM0 yr-i or m ~ 0.03, a = 0.1). Courtesy of J-M. Hure. 



7 The a-Prescription, Computation of the Vertical Structure 
and Vertical Transfer 

The computation of an accurate emission spectrum requires a model of ver- 
tical structure with the self-consistent treatment of radiative transfer. Since 
the disk is still assumed to be geometrically thin, it is again possible to decou- 
ple the vertical and the horizontal structure equations. One then solves for 
each annulus a set of differential equations in z, instead of a set of algebraic 
equations. A simplification consists, at least in a first step, of decoupling the 
computation of the emitted spectrum from that of the vertical structure, as 
the disk is optically thick and the spectrum forms in the upper layers. The 
radiative transfer is then solved using the diffusion approximation. 

7.1 Vertical Structure 

The equations for the vertical structure are: 
• the hydrostatic equilibrium, 

= -nlz ; (47) 



p dz 

• the vertical variation of the energy produced by viscous dissipation, 

^ = ^^pu{z)nl ; (48) 

• the variation of the mean optical thickness and of the surface density, 

dr , dT, 

— = —KO and — — = p , (49) 

dz dz 

where k is a frequency averaged opacity; 

• the variation of the radiation pressure in the diffusion approximation, 
assuming that the energy is transported only via radiation, 

-^^=i^(z) . (50) 
fiK dz 
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A question is raised here: how is viscous heating deposited vertically? The 
original Shakura-Sunyaev prescription was introduced to describe vertically 
averaged quantities, and it is not clear how to extend it to a z-dependant 
model. According to the a~ P prescription described above, it seems natural 
to assume that v is constant vertically, which implies that the heat deposition 
dF/ dz is proportional to p. If the opacity is dominated by Thomson scattering 
and the pressure is dominated by Prad, as is the case in the inner region of the 
disk, then Eqs. (47), ( |4^ ) and ( |50| ) lead to the result that p does not deuEiid 
on z. A different prescription is sometimes used for disks in dwarf novaeLJ. 

Other ingredients can also be added to the model, such as convective 
transport. As an illustration. Fig. |l^ displays the 2D structure of an a-disk, 
for a =|-Q.l, taking into account convective and turbulent transport and self- 
gravityEll. The parameters are those of a typical AGN (M = IO^^Mq, M = 
O.lMQ/yr, or TO ~ 0.03). Comparing the results of 2D and ID computations, 
one finds that the results are similar for the radial structure within factors of 
a few. The temperature in the equatorial plane is accurate to within SjOSo and 
the most important difference (by a factor < 4) is in the disk thicknesal3. The 
density and the surface density, the critical radius for self-gravity, the ratio 
of gas to radiation pressure ratio are correctly computed in the ID model (cf. 
Fig. H). The 2D structure avoids however some multiple solutions met in the 
ID model (cf. Fig. |). 



7.2 Disk Spectrum 

Vertically averaged a-disks, which have proven to be extremely useful in de- 
termining the radial structure of the disk, are not sufficient to determine the 
emission spectrum. This requires solving the vertical radiative transfer in a 
consistent way with the vertical structure of the disk. 

Several authors have built consistent models of the vertical structure of 
AGN accretion disks using various approximations for the-cadiative transfer. 
The first computation was performed by Laor & NetzertJ, who solved the 
radiative transfer in the Eddington approximation, with LTE, bound-free and 
free-free opacities, and with different assumptions for the stress tensor as well 
as for vertical heat deposition laws. They showed in particular that the Ly- 
man discontinuity can appear either in absorption or in emission according 
to these laws. Ross et al.E£Zl relaxed the LTE assumption andJxipk into ac- 
count Comptonization using the Kompaneets equation. StorzerE 



123 



solved the 



radiative transfer equation in 3D, also taking into account the inffuepa 
self-gravity, but keeping the LTE approximation. Shimura & Takahara 



117 



116 



also solved the vertical structure together with the transfer treated with the 
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Figure 11. Spectrum of an a-disk around a Schwarzschild BH. Solid lines represent the 
contribution from different annuli and the total spectrum (with Componization). Dotted 
lines represent the local BB approximation used for the outermost annuli. Dashed lines 
represent the total spectrum with the local BB computation. From Ross et al.^"* 



Kompaneets equation, but for a stress tensor proportional to the gas pres! 



124 



and Shields & Colemantl^ 



and not to the total pressure. Storzer et al 
addressed the question oLhow departures from LTE influence the Lyman dis- 
continuity. Dorrer et alHa solved the vertical structure self-consistently, with 
LTE and no bound-free opacities, and a pure-ILatmosphere, but taking into 



account Compton scattering. Sincell & KroliUHS solved the coupled vertical 
structure and radiation transfer with a two stream Eddington approximation, 

restricting the elements to He and H in LTE. | . 

Fig. |l^ illustrates the influence of ComptonizatioEllSJ. It shows the spec- 
tra of two (non-rotating) BHs, with different masses and accretion rates. This 
figure illustrates two interesting results: 

1. The soft X-ray excess cannot be reproduced. Although a non- negligible 
fraction of the disk luminosity is radiated in the soft X-ray range as a 
consequence of Coniptonization, it is obvious that the spectral energy 
distribution of AGN cannot be accommodated if one only assumes an 
extra power-law component exists in the soft X-ray range. Maximally 
rotating BHs produce continua which are harder those of non-rotating 
BHs, but they do not solve the problem either. 

2. The local BB gives the correct result at optical to near-UV wavelengths 
(this would not be true for the UV band and for higher massesEJ), and 
conflrms our previous conclusion that the time-delay observed between the 
optical and UV bands is too small compared with the smallest time-scale 
in the disk. 

Hubeny and coUaborator^^'ElJli solved self-consistently the vertical 
structure with an exact radiative transfer with no diffusion and escape prob- 
ability approximation, in the non LTE case, taking into account the general 
relativistic transfer function. However, they did not take into account Comp- 
tonization, or heavy elements, so their models probably do not correctly de- 
scribe the EUV. They computed the most extensive set of models to date: 
a grid of non-LTE disks around maximally-rotating Kerr and Schwarzschild 
BHs, for a wide range of BH masses and accretion rates. Some of their results 
are shown in Figs. |l2|, and 
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Figure 12. Spectral energy distribution of an a-disk around a maximally rotating BH 
(solid lines) and a Schwarzschild BH (dashed lines), with M = 10® Mq and various values 
of the accretion rates. The curves correspond to -L/-LEdd= li 1/2, and 1/4, etc., for the 
Kerr BH; and the corresponding values for the Schwarzschild BH are 5.613 times larger. 
Prom Hubeny et al.*^ 



Figure 13. Optical/ultraviolet spectral slope between 1450 and 5050A for an a-disk 
around a Kerr BH with 37° viewing angle, and a = 0.01. Solid curves: computed models; 
dashed curves: local BB. From top to bottom, the curves are computed for 9 values of M: 
1/8, 1/4, 1/2, 1, 2, 4, 8, 16, and 32 xIO^Mq. From Hubeny et al.*^ 



Figure 14. Lyman edge strength as a function of Eddington ratio for an a-disk around a 
Kerr BH with M = 10^ Mq . The curves are labelled with the value of cosi. From Hubeny 
et al.« 



Fig. compares the energy distributions of disks of Schwarzschild 
and Kerr BHs of 10^ Mq viewed at 37°, for different disk luminosities. 
Schwarzschild and Kerr BHs produce quite similar optical spectra, except 
that in Kerr BHs the spectrum is shifted towards higher frequencies. 

Fig. ^ displays 1450 to 5050 A spectral slopes, (3 (defined here as F^, oc 
ly^), as a function of the Eddington ratio for Kerr BHs of different masses. 
The local BB models depart from real transfer models except for low accretion 
rates. This is due to the relatively high masses considered. Bath sets of models 
can reproduce the mean observed colors of the LBQS samplec3. However, one 
should note that the models cannot account for the positive slopes observed in 
the majority of Seyfert galaxies. 

Finally, Fig. [l^ shows the strength of the Lyman edge in a Kerr BH 
with different disk inclinations. The Lyman discontinuity is always more 
intense than observed, except for very low accretion rates or very large disk 
inclinations, not plausible in the framework of Unification. It is however 
possible that the addition of metal opacity (not taken into account in the 
computations) changes this result. 

This discussion illustrates the fact that the computation of disk spectra is 
difficult and involves poorly known parameters. Besides, several other param- 
eters not mentioned here complicate the problem, such as the value of a, the 
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Figure 15. Sketches of several models of irradiated disks (the lamppost, the disk-corona, 
the patchy corona), and of the cloud model. 



inner and outer radii of the disk, the different assumptions concerning heat 
deposition, etc. Therefore, in some cases it may not be a bad choice to use 
the simplest and least model-dependent solutions. Despite these uncertainties, 
the following results can be considered as secure: 

• the spectral energy distribution does not fit the observations in the soft 
X-ray range: another emission mechanism is required; 

• the disk spectrum in the optical to near-UV range is too flat compared 
to that of many Seyfert nuclei; 

• the computed Lyman discontinuity is too strong or has to be fine-tuned; 

• the small or absent time-lag between the optical and UV bands is not 
compatible with the viscous time-scale, or even with the dynamical time; 

• the disk is not emitting hard X-rays at all. 

Additionally, we paid no attention to polarization properties. They are 
also hard to explain in the context of the standard disk model. 

8 Other Disk Models 

8.1 The "Lamppost" Model 

The discussion on the standard disk model has shown that the regions emitting 
at UV and optical wavelengths are located at large distances from each other, 
so the small time-lags between the (well correlated) optical and UV light 
curves require a causal link propagating at the speed gfJight. This led to 
questioning the validity of the thin steady disk modelEi'LiI, and to propose 
that the disk was actually irradiated by the X-ray continuum, j-aiid jemitting 
as a result of both gravitational and externaJ-radiative heatingEZlH'E£3. The 
discovery of the iron line and the X-ray humpO implied reprocessing by a cold 
medium, in agreement with this model. 

There then followed the era of "irradiated disks" . These models included 
an X-ray point source located at a given height A above the center of the disk, 
or a spherical source of radius A around the BH. The former is now referred 
to as the "lamppost model" (cf. Fig. ^5|). 

The flux falling on each face of the disk i^x can be expressed as 

where /x is the fraction of the bolometric luminosity emitted by this source, 
6 is the angle of the light rays to the normal, and a — \/2 for a point source 
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and a — 2/37r for a spherical source. The second approximation is vahd for 
i? ^ A, i.e. for the optical to near-UV emitting regions, located much further 
away than the X-ray source. X-ray variability studies suggest that A ^ 10i?G- 
Therefore, the external flux varies with the radius like the viscous flux. This 
external irradiation has several effects on the disk: 

• A major fraction of X-rays is absorbed by photoelectric processes in the 
upper layers of the atmosphere, and this creates a hot and ionized "skin" 
which emits the iron line. 

• The remaining fraction is Compton scattered in this ionized atmosphere 
and gives rise to the X-ray "hump" . 

• The absorbed fraction is reprocessed into UV, of which about 50% is re- 
emitted outwards and 50% towards the interior. This latter fraction can 
modify the internal structure of the disk. The effect is however negligible, 
unless /x ~ 1- 

Indeed, the condition for radiative heatiug to dominate the viscous heat- 
ing when the disk is optically thick isoa 

Fvisc (1 + ^tr) < ^ . (52) 

Since this condition is not realized, the radial structure is not modified 
(unlike in the disk-corona model). 

• The fraction of UV radiation re-emitted outwards is added to the emission 
due to the viscous flux. 

The emitted spectrum has been copaputed using different schematiza- 
tions. Except for the most recent onea^a, the models consist of a disk of 
constant density in the vertical direction, irradiated on the top by a power- 
law continuum and on the bottom by a BB (or by a Wien spectrum) to 
approximate the cfEciitpOf jnjsmus heating. The transfer is treated with vari- 
ous approximationgl2aJl3x3'E£il or with a Monte Carlo method coupled with 
a photoionization codaHj. An example is given in Figs. |l|and|l|. Fig. |l6| 
shows that iron is strongly ionized in the upper layers. The spectrum dis- 
played in Fig. OJJlustrates the effects mentioned above. It is computed with 
the code TITANEZI, which solves the complete non-LTE transfer inside lines 
and continuum, with photoionization and thermal equilibrium, and is coupled 
to a Monte Carlo code to handle Compton diffusion. The strong UV-bump 
corresponds to the viscous release, to which a small (~ 10%) contribution 
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Figure 16. Ionization state of iron in the upper layers of an irradiated disk. The model is 
a slab of constant density lO^'^cm"^, irradiated on one side by a power-law Fy oc u~^'^ 
from leV to 100 keV and on the other side by a BB with T = 10^ K to mimic the viscous 
heating. The ionization parameter is ^ = 300, and the luminosity of the power-law is equal 
to the luminosity of the BB. Computations made with TITAN. 



Figure 17. Emission spectrum for the same model as Fig. |l6| The dotted line represents the 
underlying transmitted BB spectrum; the dashed line, the reflected spectrum; and the solid 
line, the total spectrum emitted by the disk. The spectral resolution (R) of the spectrum 
has been degraded to 30. 



of reflected radiation is added. It is interesting to note that the Lyman dis- 
continuity, present here in absorption, is quite weak, though not completely 
absent, as in the observed spectra. In the soft X-ray range the spectrum has 
a slope of the order of 1.5 in flux, in agreement with the observations, and 
displays many lines and edges from a highly ionized medium. SujcIi lines are 
observed in the Chandra spectra of NGC 3783^, and NGC 5548y, unfortu- 
nately strongly blended with absorption lines from the Warm Absorber, and 
perhaps partly produced by it. This continuum is the "reflected" spectrum, 
though obviously it is not due to pure reflection! In the hard X-ray range 
an intense iron line is produced at 6.4 keV, due to Fe XVII and less ionized 
species. 

This model is able to solve several problems raised by the standard non- 
irradiated disk model, like the presence of the soft X-ray excess, the iron line 
and of the "hump" ; but it retains a few unsolved problems: 

• Only a small fraction of the UV flux can be produced in this way as the 
X-ray luminosity is smaller than L(BBB). One deduces that the optical 
to UV emission is only slightly modiflcd with respect to the standard 
disk (see also Fig. |l^). Since the local BB picture holds for the optical 
to near-UV emission, one has Tbb = [{Fyisc + ^x)/cr]^/'*, and therefore 
Tbb oc except in the very central regions, where anyway the BB 
approximation is not valid and only EUV is emitted. Consequently to 
get a steep spectrum, F^, cx i^"" with a > as observed in Seyfert nuclei, 
one should seek another geometry for the irradiation of disks. 

• The disk will respond to a variation of the central source on a time-scale 
of- 1.2\*J^rr^/^f{Ry^'^m^/^A4^'^ days, according to Eq. (|2l|). Thus, 
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rapid X-ray fluctuations will be erased by the light travel-time towards 
the disk, as observed. But the UV-optical light curves should always 
lag behind the X-ray one, while the opposite is sometimes observed, for 
instance in the case of NGC 7469. Also the model cannot explain the 
fact that the UV light curve has the same variability amplitude as the 
X-ray with the characteristic UV time-scale. 

• The problem of the Lyman discontinuity remains, although improved 
with respect to non-irradiated disk. 

• Finally, the model requires the presence of an ad-hoc X-ray source for 
which no physical emission mechanism is proposed. 



8.2 What Can We Learn about the Emitting Medium from the X-ray 

Spectrum?: Where Compton Scattering Makes its Appearence Again 

We know that the X-ray spectrum is made of a reflected part (the "cold" 
iron line, the hump), and of a primary spectrum. The reflection takes place 
on a "cold" medium, which most probably is that emitting the UV-bump. 
The correlated variations and the small time-lag of the UV and X-ray fluxes 
confirm that the two emitting media are "seeing" each other. The high energy 
cut-off of the primary spectrum at about 100 keV tells us that it is emitted by 
a hot medium with an electron temperature such that — kTe/mc^ 0.1 
(if the emission is mainly due to electrons, which is most likely). 

Three thermal mechanisms can compete in this hot medium: free-free, 
cyclotron, and inverse Compton emission. The time-scale for inverse Compton 
cooling is smaller than the time-scale for free-free cooling when the following 
inequality holds 

^ ^ Lw ^ ^ Lw 1 ^ 0.039,-1/2 ^ 10-2 ^ (53) 
R mc'^ L^dd r 

where I is the "compacity" of the UV emission region. 

The dimension of the UV emission region is typically 100_Rg. One deduces 
that the condition is easily fulfilled, so inverse Compton cooling dominates 
over free-free cooling. It dominates also over cyclotron cooling unless the 
magnetic pressure is of the order of the density of radiation, which implies 
a very high value of the magnetic field. In conclusion: the optical- UV and 
the X-ray emissions are probably strongly coupled through direct (in the cold 
medium) and inverse (in the hot medium) Compton scattering. 

In the hot Thomson thick medium there were multiple scatterings of hard 
photons. In the hot Thomson thin medium with soft photon input there is-a 
small number of scatterings. This model was proposed for X-ray binariestij. 
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The mean amplification per scattering is A = 1 + 48e — hi'o /rrieC^ . If tt < 
1, the probability of m scatterings is r™. After m scatters, v/vq = A™. This 
creates a series of humps with decreasing intensities. One can thus approxi- 
mate the Comptonized spectrum as a power- law /(i^) — I {vq){v / v^f'^'^'^'^^ / \nA_ 
For a given energy of the soft photons, the spectral index is therefore defined 
by Tt and T^. Moreover no photons can be scattered to an energy greater 
than that of the thermal electrons, and thus, one expects the power-law to be 
cut at hf/kTe ~ 1. 

Since the cooling of the hot gas is due to inverse Compton scattering, its 
temperature can be determined from an energy balance equation, provided 
that the heating is known, and thus one can completely determine the spectral 
distribution. This led to the introduction of a disk model where the release 
of gravitational energy occurs in a hot optically thipL-'Japmna" above the disk, 
through buoyancy and magnetic field reconnectionl£a'E£l'LZl. If a large fraction 
(/ ~ 1) of the gravitational release takes place in the corona, the disk is almost 
exclusively heated by the backscattered radiation coming from the corona. 
It reprocesses this radiation into a soft BB, and the corona upscatters this 
radiation while being Compton cooled by this process (cf. Fig. |l^). This 
radiative coupling between the disk and the corona and the energy balance of 
the corona impose a relation between Te and tt (the smaller tt, the larger 
Te), such that the X-ray spectral index is close to unity. Since an increase in 
Tt causes both a decrease of Te and a steepening of the spectrum, a relation 
between the cut-off at high energy and the spectral index in the 2-10 keV 
range is predicted (cf. Fig. [Ts] ). 

Since the ratio between the inverse Compton luminosity and the soft lu- 
minosity is of the order of unity, the model predicts a BBB/X-ray luminosity 
ratio smaller than observed. This is why a variant, of the disk-corona model, 
the "patchy corona" (cf. Fig. 15), was introducedEa. The corona is not homo- 
geneous, but is made of a few blobs, which could be due to the formation of 
magnetic loops storing energy and releasing it rapidly through reconnection 
like in solar flares. Below a flare, the disk is heated to a higher tempera- 
ture. The rest of the disk radiates through viscous processes. An important 
difference with the continuous corona is that the fraction of the disk located 
below the "flare" is much more irradiated and consequently its upper layers 
are strongly modified, in particular due to the large radiative pressure. 



8.3 The Two-Phase Disk 

In the previous model the gravitational release takes place mainly in the 
corona. This implies strong structural changes of the disk, both vertically 
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Figure 18. Comptonized spectra for different couples of (tTiSg) [(0.63, 0.11), (0.32, 0.21), 
(0.2, 0.3), (0.1, 0.5)]. The soft photons correspond to a BB with /cTbb =100 eV. From 
Haardt et al.^^ 



and radially, with respect to the standard model 

Disk- 
papersEll' 



ina 1 1 models. have been studied in many 

3'liH3'E£3'E2j'lH2l. Since the disk has no or only weak inter- 
nal heat generation, it is almost isothermal and has a weak radiation 
pressure. It is therefore very dense and supported only by gas pressure. 
Consequently it is not subject to the viscous/thermal instability like radiation 
supported disks. This is a great advantage, but on the other hand, since 
these disks are very dense, they are gravitationally unstable at all radii for 
relatively modest luminosities. Other differences with the lamppost model 
is that the corona pressurizes the disk, and there is conductive transport 
between the disk and the corona. 

Generally, the corona is a two temperature plasma, where ions are heated 
and transfer their energy to the electrons through Coulomb collisions, and the 
electrons are Compton cooled. The corresponding equilibrium temperature is 
close to virial for the ions (lO^^K), and of the order of lO^K for the electrons, 
as already mentioned. Therefore, the gas pressure supported corona is at the 
limit of being geometrically thick, and sometimes it is quasi-spherical. Some 
people think that a collective process able to rapidly transfer energy between 
ions and electrons can maintain the two constituents at the same temperature. 

In these models, the fraction of energy generated in the corona is a free 
parameter given for instanceES by a power-law in R. By changing this fraction 
and the geometry of the hot region, the models can fit the observed X-ray 
spectra. But one should search for a physical basis for the choice. 

Some models take into account mechanical transport only in the radial 
direction, and others add the vertical transport by conduction or evaporation 
between the disk and the corona. Models with the local corona strength, cal- 
culated on the basis of the disk evaporation rate, tend to predict a relatively 
weak rnrnna ahnvp a disk, with possible total disk disruption in its inner- 
most partllij'lll These results are supported by 2D numerical time-dependent 
computationsO. However, these models do not predict strong the soft X-ray 



excesses present in composite spectraEj, as Comptonization in the corona is 
too weak and the emission from the innermost single phase flow arises in the 
hard X-ray part of the spectrum. 
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8-4 The Iron Line: a Proof for Irradiated Cold Relativistic Disks? 

When an X-ray photon is absorbed to eject one of the two K-shell electrons 
of an iron ion, this photoelectric process is followed by the de-excitation of an 
electron from the L-shell, and then, either by the emission of a "fluorescent" 
K-photon, or by a subsequent Auger ionization (escape of an electron from 
an upper shell without emission of radiation). The probability of these two 
processes is of the same order and depends slightly on the ionization state of 
the ion. The energy of the fluorescent photon is almost constant and equal to 
6.4 keV for Fe I to XVII, then it increases to 6.9 keV for the H-like Fe XXVI. 
Actually this transition and the one corresponding to the He-like Fe XXV are 
not "fluorescent" but simply normal resonance La transitions. Besides, from 
Li-like (Fe XXIV) to H-like iron, the Auger effect cannot occur, owing to the 
lack of two electrons in the L-shell. The absorption threshold is 7.1 keV for 
photoionization of a neutral atom, and it increases to 8 keV for Fe XVII, to 
reach 9.3 keV for Fe XXV. 

The ionization state of iroii |d£pends on the ionization parameter ^. For 
constant gas density one findsO^O: 

for ^ < 100: a "cold" line is produced at 6.4 keV by ions less ionized than 
Fe XVII, and the iron K-shell absorption edge is small since the contri- 
bution of the other elements to the photoelectric opacity is large. 

for 100 < ^ < 500: iron is in the states Fe XVII to Fe XXIII. There is a va- 
cancy is the L-shell of the ion, so Ka photons are resonantly absorbed 
and trapped until they disappear by the Auger effect. Only a few line 
photons can escape and the Fe K line is very weak. The lower-Z elements 
are ionized, leading to a moderate iron absorption edge. 

for 500 < ^ < 5000: the line photons are still subject to resonant scattering, 
but since they are not destroyed by the Auger effect they can escape the 
disk and produce an intense "hot" iron line at 6.8 keV. The absorption 
edge is strong. 

for ^ > 5000: iron and the other elements are completely ionized, so there is 
no line and no edge. 

As we shall see in the next section, this is no longer valid for a constant 
pressure medium or for a disk in hydrostatic equilibrium. 

The flux of the line and the strength of the ionization edge depend on 
the ionization state, but also on the iron abundance, the geometry of the 
reprocessing medium, and, in the case of a plane-parallel slab, on the incli- 
nation of the line of sight and the illuminating source itself. We shall not 
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Figure 19. Profile of the iron fine emitted by a disk around a Schwarzschild and a maximally 
rotating Kerr BH. The disk is assumed to extend down to the last stable orbit. From Fabian 
et al.'^^, ©2000, Astr. Soc. of the Pacific, reproduced with permission of the editors. 



give here more details about these dependences (see Fabian et al.0). The 
observed mean intensity, centroid energy, and equivalent width, lead to the 
conclusion that the line is produced in a relatively "cold" medium (meaning 
that iron is less ionized than Fe XVII), covering about 27r steradians of the 
central source (therefore an infinite slab irradiated by an X-ray source is an 
acceptable geometry). 

It is widely accepted that the line profile is a proof of the formation of the 
line at the surface of the accretion disk, very close to the BH. The extended red 
wing is then due to gravitational redshift, while the blue wing is relativisticaly 
boosted due to the large rotational motion. As an illustration. Fig. [T9| displays 
on the left the profile of the line emitted by a disk around a Schwarzschild BH 
extending from 6 to 30i?Gi as a function of the inclination of the disk. The 
intensity of the blue wing depends strongly on the inclination. The figure on 
the right compares the profile of a line formed in an accretion disk around a 
Schwarzschild and a Kerr BH, for an inclination of 30°. In the case of the 
Kerr BH, the disk is assumed to extend down to 1.25i?Gi so the gravitational 
redshift is stronger and the red wing is more extended. 

Fitting the profile of the line in several objects led to a number of surpris- 
ing results. First the absence (or the weakness) of the blue wing corresponds 
to an inclination close to 30°. Second, at least in the two best observed ob- 
jects, NGC 3615 and MCG-6-30-15, the disk should extend down to 2Rg to 
account for ±hn broad red wings. Fabian et al.Ell argued against Reynolds 



& Begelmar£23 for a Kerr geometry, as the gas falling inside the last stable 
orbit on a Schwarzschild BH would be highly ionized and the spectrum would 
display a strong absorption edge, inconsistent with the observed data. 



. 5 More Realistic Computations of the Emission Spectrum 
A first attempt to compute the emission of a_disk topped by a corona in a 



consistent way was made by Sincell & KrolikE£ll. Their transfer method was 



already mentioned in the discussion of the standard disk. They relaxed the 
condition of constant density, and solved the coupled vertical structure (i.e. 
hydrostatic equilibrium of the irradiated layer) and the radiation transfer, 
assuming that the energy release takes place entirely in the corona. The 



suzy-collin-mexico: submitted to World Scientific on February 1, 2008 30 



Figure 20. Comparison between the temperature structure of an irradiated layer for con- 
stant density and for constant pressure. The curves are labelled with the log of the flux 
incident on the layer. The constant density layer and the surface of the constant pressure 
layer have a density lO'^^ cm""^, the incident spectrum is a power-law with a spectral index 
unity, from 0.1 eV to 100 keV. Computations made with TITAN. 



density in the X-ray heated skin is very high (~ 10^'* cm^'^) and justifies their 
LTE approach. Integrating the emission over the radius, they computed the 
UV continuum emitted by the X-ray heated skin. Their calculations predict 
a strong Lyman discontinuity in emission, in strong disagreement with the 
observations . _ _ 

Nayakshin and coUaboratorsE^'Ej also solved consistently the radiation 
transfer, with a complete treatment of the ionization and thermal equilibrium 
(not limited to H and He) and the hydrostatic equilibrium of an irradiated 
disk, not restricted to the case of a strongly pressurizing corona. They stressed 
that hydrostatic equilibrium is extremely important for the structure of the ion- 
ized skin, owing to the the thermal instability mechanism linked with the 5 
vs. T curve (cf. Fig. ^): when the temperature reaches the unstable part 
of the curve, it "jumps" to the other stable part. In contrast with the con- 
stant density case, which predicts a smooth decrease of the temperature with 
increasing depth inside the irradiated layer, hydrostatic equilibrium predicts 
a sharp transition between the hot and the cold stable phases. The same 
phenomenon occurs in the constant pressure case, as shown in Fig. One 
also sees that the hot layer is thicker in the constant density case, owing to 
the rapid increase of density. As a consequence, for constant pressure the 
spectrum displays a UV-bump restricted to a one temperature BB, while for 
constant density the UV-bump is much more extended in frequency, and it 
joins smoothly the soft X-ray range (cf. Fig. 21). Fig. shows that for con- 
stant density, no "cold" iron is present, unlike in the constant pressure case, 
where the cold layers are closer to the surface, so the soft photons emitted by 
these layers are able to leave the medium. 

For hydrostatic equilibrium, the cmissivity of the line and the spectral 
distribution of the continuum reprocessed in the atmosphere of an irradiated 
disk depends strongly on a dimensionless "gravitational parameter" A r 
portional to the ratio of the gravity to the radiation pressure tcrmE^'M 
Fig. ^ shows that for constant density, when the iron line is intense, it is 
always highly ionized. For hydrostatic equilibrium (which isi referred to as 
the self-consistent case), when the iron line is intense, it is a "cold" line. This 



suzy-collin-mexico: submitted to World Scientific on February 1, 2008 31 



Figure 21. Comparison between the spectrum reflected by an irradiated layer for constant 
density and_for constant pressure. Same density and incident continuum (the dotted line) 
as in Fig. 



Figure 22. Zoom of Fig. U^. Computation made with TITAN and NOAR. 



Figure 23. Comparison of the constant density case with the hydrostatic equilibrium case 
("self-consistent case"). The curves are labelled with ^ in the first case and with the 
gravitational parameter A in the second case. One sees that the iron line is "cold" in the 
self-consistent case, and "hot" in the constant density case. From Nayakshin et al.^^ 



might explain the predominance of "cold" iron in AGN spectra. The observed 
correlations are better accounted for by a strong illumination of the disk, cor- 
responding to the "patchjj" or "flaring" corona model, than by the lamppost 
or the continuous coronaO. Clearly more of these self-consistent computations 
(unfortunately enormously time consuming) are needed to get a real physical 
understanding of the complex observed behavior of AGN spectra. 

9 An Alternative to the "Cold Relativistic Disk"? The "Cloud 
Model" 

The iron line and UV continuum are both due to reprocessing of the X-ray 
continuum, for instance in the atmosphere of an irradiated disk. Hence the 
regions contributing the most to the variable fraction of the UV flux should 
give rise to an iron line, and inversely the region emitting the iron line should 
produce a UV or a EUV excess of continuum emission. 

We have seen that the optical- UV variable continuum is emitted at large 
distances from the BH (at B ~ 100 to 1000 Rq), with the result that a fraction 
of the iron line is also reprocessed in this region. The iron emission will appear 
as a narrow line at 6.4 keV, constant on UV variation time-scales. It is not 
clear whether a narrow core is always observed, but we know that the broad- 
line stays constant on time-scales of the order of days while the underlying 
X-ray continuum varies. Also there seems to be a contradiction between the 
short time-lag of the optical and UV light curves and the relative large time- 



suzy-collin-mexico: submitted to World Scientific on February 1, 2008 32 



lag between the UV and the X-ray hght curves. Generally speaking both 
the iron line and the optical-UV-X observations, in particular their variability 
properties, have been surprising in the context aLthc-disk model. 

An alternative geometry has been proposedE3'llj'l23, where the Fe K line 
and the UV continuum may arise from dense clouds (optically thick for free- 
free absorption, but not Thomson thick) with a large covering factor embedded 
in a hot quasi spherical medium, rather than from a disli_.A variant of this 
model, with Thomson thick clouds, has also been proposedlia. Fig. displays 
a sketch of the "cloud model". In this model Compton scattering broadens 
the iron line. 

A physical basis for the model can be found in a new inhomogeneous 
equilibrium of the disk where the bulk ofJjhe mass rests in dense clumps 
magnetically connected in a bath of hot gaa£3. One advantage of the model is 
that it is thermally and viscously stable. The transport of angular momentum 
is insured by collisions between clumps. The requirements for dynamical and 
thermal equilibrium determines the parameters of the system, for example the 
temperature of the hot gas and the optical thickness (quite similar to that of 
the corona in the disk-corona model) of the clumps which are Compton thick 
with a covering factor close to unity. 

Another possibility exists in the framework of a nuclear wind. There 
are several theoretical arguments for the existence of hydrodynamic winds or 
radiation-driven flows close to the BH. On the other hand, the fact that the UV 
and soft X-ray absorption lines are blueshifted is observational evidence for an 
outflowing medium. One can thus imagine that the quasi-spherical medium 
invoked in this model is simply the basis of the Wapn Absorber, for instance 
similar to the wind envisioned by Murray & ChiangE^. One interesting aspect 
of such a model is that the iron line would be easily redshifted by the Doppler 
effect from the position of the Fe XXV or Fe XXVI hne (6.7 keV) to 6.4 
keV or less, since the observer would be seeing a line emitted only from the 
illuminated side of the clouds. ___ 

This model is able to account for the UV-X spectrum of AGNta'Ea'N . The 
observed X-ray continuum is then a mixture of the primary continuum and of 
the continuum reflected by the illuminated side of the clouds, while the BBB is 
produced both by intrinsic emission of the clouds heated by the X-ray source 
and by reflection on the illuminated side of the clouds. Similar to the disk- 
corona model, the X-ray source can be the result of Compton upscatterings 
of the UV photons emitted by the surrounding clouds. The radius of the cold 
medium is r ~ 40mn^^^^ ^^^^3 where ni4 is expressed in units of 10^^ 
cm~^ and ^3 in 10'^. Fig. ^ shows an example of the spectrum obtained for 
such a model. 
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Figure 24. Spectrum emitted by a spherical shell of clouds surrounding a hot medium that 
upscatters the UV photons emitted by the irradiated clouds. The covering factor of the 
clouds is 0.9, and the ratio of the radius of the hot and cold media is 0.4. Spectrum 
computed with TITAN and NOAR. The noise in the right-hand side is a numerical artifact 
of the Monte Carlo method. The spectral resolution of the left-hand side is 100. 



Figure 25. Amplification factor due to multiple reflections in the cloud model, assuming 
that the X-ray source is point-like. Starting from ^ = 3000, the effective primary ionization 
parameter illuminating the clouds reaches § = 10880. The covering factor of the clouds is 
0.9, and the primary spectrum is a power-law with a spectral index equal to unity between 
0.1 eV and 100 keV. From Abrassart & Dumont'^. 



Figure 26. Fit of the mean iron line profile with the cloud model. The incident spectrum 
on the reprocessing shell is a power-law with unity spectral-index between 0.1 eV and 100 
keV, the covering factor is 0.9, and the shell contains highly ionized (t; = 10000) and nearly 
neutral = 300) clouds in equal proportion. From Abrassart &; Dumont'^. 



In this closed geometry the UV photons which are not upscattered to the 
X-ray range by clouds in the hot medium, undergo multiple reflections which 
amplify the UV-bump. Consequently the clouds "see" a spectrum richer in 



soft photons than the observer (cf. Fig. 25). This leads to a small Comp- 
ton temperature of the clouds, which is important to account for Compton 
downscattering of the iron line photons (cf. below). |— , 

Such a medium has the power to Compton-broaden the iron lineE^. There 
are actually two possibilities^: 

• the hot medium is surrounded by a cold Compton thick medium, so the 
iron line and the hump are eniitted by the Comptonizing medium itself 
(Comptonization by reflection^; 

• the iron line is intrinsically narrow because it is emitted at ^ 100i?G 
from the BH, and it is Compton broadened when crossing a surroppding 
ionized Compton thick medium (Comptonization by transmission )Eil'E3u. 

An example of the model fit of the mean AGN spectrumil with Comp- 
tonization by reflection is shown in Fig. |2^. Two types of clouds arc required: 
"cold" clouds, which emit the "cold" component of the iron line, and highly 
ionized clouds, which injure Comptonization of the iron line. The BBB is 
produced by both mediaEI. 

An interesting aspect of Comptonization by transmission is that it can 
account not only for the broad red wing of the iron line, but also for the 
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Figure 27. Profile of an intrinsic cold narrow 6.4 keV iron line, transmitted through an 
extended shell of Thomson depth 3, with a density varying as and § = 2 X lO". The 

incident spectrum is a composite AGN^^. From Abrassart & Dumont*. 



Figure 28. Same as in Fig. |27|, but showing the effect on an extrinsic Lyman edge in 
absorption. From Abrassart 



absence of a Lyman discontinuity in the UV spectrum without appeahng to 
a fine tuning of the parameters (cf. Fig. 



The transmission model has been rejecteclEI'E21l based on the size given 
by variabihty. The very high ionization state required to avoid an intense 
iron absorption edge imphes too large a dimension, in particular in the case 
of MCG-6-30-15. However, assuming that the Comptonizing medium is ho- 
mogeneous (the most conservative case), one deduces from Eq. (|^) assuming 
^ = 10® and tt ~ 3 a dimension of 0.03 It-d, consistent with variability time- 
scales. It corresponds to a distance of about 100i?G if the BH in MCG-6-30-15 
has a mass of a few IO^Mq, as it is now believed (note that such a small mass 
implies that this object is radiating close to its Eddington luminosity, making 
it peculiar among Seyfert nuclei). However a serious argument against the 
transmission model, is that the energy cut should take place at 500 keV / t^, 
i.e. at about 50 keV, while it seems to be observed above 100 keV. 

A strong argument against the reflection model is the fact that it still has 
not been shown able to account for a very extended red wing with no blue 
wing, as observed in MCG-6-30-15. However the parameter space is far from 
having been completely explored, as these computations are extremely time 
consuming. 

The main difference with the disk model is that the region where the iron 
line is emitted has not to be located close to the BH. This is in agreement 
with the absence of short term variations of the line strength correlated with 
variations of the continuum (though the rapid variations of the profile are 
difRcult to account for except with an ad hoc geometry). It allows also for the 
UV and optical emission to arise at similar distances, of the order of 100i?G- 
On the other hand, in this model the clouds can randomly obscure the line of 
sight, thus accounting for short term variations of the primary X-ray sourceu. 
Another major difference with the disk model is that it also explains naturally 
a large Luv / Lx ratio, owing to the large covering factor of the "cold" medium. 
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Note that this model is in a sense "non local" , as the cloud system can be 
compared to a photosphere which is reprocessing the gravitational energy that 
is emitted much closer to the BH. |_. 

Finally an intermediate solution has been proposedE^ where a system of 
quasi-spherically randomly keplerian moving clouds is located relatively close 
to the BH. In this scenario the iron line is broadened both by relativistic 
effects and by Comptonization. 

10 And What Happens at the Largest and at the Smallest 
Accretion Rates? The Case for Thick Disks. 

Up to now we have concentrated on accretion rates similar to those observed 
in Seyfert nuclei and modest quasars, where geometrically thin disks are a nat- 
ural solution of the a-prescription. It is interesting to consider what happens 
at both higher and lower accretion rates, where thick disks are predicted. 

In thin disks, vertical motions and radial pressure gradients are negligible, 
and it is therefore possible to neglect the flux advected towards the center. 



When the disk is geometrically thick, the radial velocity is comparable to the 
rotational velocity, and this induces strong radial pressure and temperature 
gradients. The advected flux must then be taken into account. 



where i^cooi is the heat flux radiated away. Let us consider the two extreme 
cases: 

I. for m or ^ I: Owing to the large accretion rate, the flow is dense 
and optically thick, the inner regions are thus supported by the strong 
radiation pressure. Consequently, the disk is inflated and becomes ge- 
ometrically thick. Heat advected to the BH exceeds the heat radiated 
away, and 77 decreases with m □. For super-Eddington accretion rates 
these disks have relatively small luminosities, of the order of iEdd- They 
radiate mainly in the soft X-ray range and it is suggested that they could 
account for Narrow Line Sy Is, which probably pcadiate near their Ed- 



2. for m < a^, i.e. typically < O.OI: Owing to the small accretion rate, 
the flow has a low density and does not radiate efficiently. Like the 
coronae discussed previously, the protons remain very hot and the disk 
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Figure 29. Spectrum of an ADAF with a cold disk. The 4 curves correspond to different 
positions of the transition between the ADAF and the cold disk: infinity (no cold disk, pure 
flat ADAF spectrum), fOOO Rq, IOO-Rq, IO-Rq. The importance of the blue bump and its 
peak frequency increase when the transition radius decreases. Courtesy J. -P. Lasota. 



is supported by gas pressure, becoming geometrically thick and optically 
thin. Several cooling processes can dominate wtpJJwf-Qr not soft photons 
are present to induce inverse Compton cooling^ElO'OL^. Since these disks 
are not efficient radiators, a large fraction of the gravitational energy 
will be advected towards the BH and simply "swallowed" . These disks 
are referred to as "Advection Dominated Accretion Flows" (ADAFs). 
Both solutions (i.e. the standard disk and the ADAF) exist, without 
any understanding on how to switch between them. If the ADAF is 
surrounded by a cold disk the emission spectrum will depend on the 
location of the transition between the cold disk and the ADAF (cf. Fig. 

). In this case the hot medium constituting the ADAF can be cooled 
by inverse Compton scattering of the soft photons produced by the cold 
disk. If a nuclear radio source is present, the synchrotron photons can 
play the same role in cooling the hot gas. 

Much recent attention has been devoted these last years to ADAFs. They 
are generally thought to be present in low luminosity objects, such as LIjN,- 
ERs and the nuclei of elliatical galaxies (M ST^B), the nucleus of NGC 4258Ea, 
and the Galactic CenteiO. All show a non-resolved X-ray source but no opti- 
cal/UV continuum, as one expects from an ADAF. The theoretical possibility 
of an ADAF is a matter of curijent debate, and observationally it is contro- 
versed (in particular for Sgr A*Q). 

11 A "No Man's Land": the Disk between 10^ and 10^ 

We have seen that the disks become gravitationally unstable beyond a few 
10^i?G- Moreover, the viscous time-scale given by the a-prescription for trans- 
porting the gas towards the BH at such distances is larger than the lifetime of 
the active nucleus. We know that further from the center, at about 10^ Rq, the 
disk is globally gravitationally unstable and the supply of gas can be achieved 
by gravitational torques (cf. F. Combes in these lecture notes). However we 
know nothing about the structure of the disk in the intermediate region and 
how the outward transport of angular momentum is ensured. 
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Several solutions have been proposed. A magnetic field anchored in the 
disk can transport angular momentum non-locally, and possibly prevent grav- 
itational instability. Another solution consists of a disk made from marginally 
unstable and randomly moving clouds with large bulk velocities, where the 
transfer of angular momentum is provided by cloud coUisionsO. A more elab- 
orate model is proposed for NGC 1068l£2l in which a disk of gas clumps under- 
goes mutual gravitational interactions. These disks would be geometrically 
thick, and are actually similar to the tori invoked in the Unified Scheme. 

The observations of purely keolErian velocities in masers lying at about 
10^ Rg from the BH in NGC 425a£Zl argue strongly for a geometrically thin 
warped disk. Collin & ZahnEEl have therefore adopted the view of a marginally 
stable disk where unstable clumps collapse until protostars are formed. (How- 
ever, such a picture is only well suited for low Eddington ratios. For large 
Eddington ratios the number of newly formed stars is high, and the disk will 
become highly turbulent and inhomogeneous.) Stars embedded in an accre- 
tion disk should accrete at a high rate and rapidly acquire masses of a few 
tens of M0 , leading to a "starburst" of massive stars (not comparable in size 
with real starburstsp. The mass transport can thus be driven by supernovae 
or by stellar outflows, which at the same time can contribute to the heavy 
element enrichment observed in the BLRs at high-z and in BALs. A larger 
scale starburst could also be induced in the shocked gas of supernova rem- 
nants. Several pieces of evidence have been discussed in these lectures for the 
existence of starbursta in the very central region of Seyfert galaxies. 

Finally OstrikeiH made a very interesting suggestion that has never been 
seriously considered. He stated that if the disk is embedded in a massive 
stellar system, interactions between the stars and the disk can remove angular 
momentum from the disk at the rate expected for an AGN. 

12 Conclusion 

It is not easy to conclude such a chapter because the results are really diffi- 
cult to summarize. We have seen that the tremendous increase of data with 
the advent of X-ray missions such as ASCA, RXTE, and now Chandra and 
XMM, has considerably complicated the problem, and now precludes any 
simple interpretation. We have seen that the most detailed models, used to 
make observational predictions, are still based on crude approximations, such 
as geometrically thin, stationary and axisymmetric disks. For instance the 
models are generally assumed to be steady-state, despite the evidence that 
there exist strong variations on small time-scales. The viscosity is currently 
parametrized by the a-prescription, but other prescriptions are not excluded. 
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There are heated debates about the vahie of a, whether it is constant or not, 
whether the stress tensor is proportional to the total pressure or to the gas 
pressure. We don't know the influence of the magnetic fleld. We have no idea 
of the geometry of the inner regions. We know that there is a quasi-spherical 
hot accretion flow, but at the same time we believe that the cold thin disk 
extends down to a few gravitational radii from the BH, and in this case we 
know that it should be imstable. The status of the hot flow is not clear yet. 
In short, one could say that the more sophisticated the models, the larger the 
uncertainties! To continue to make progress, we must obtain very good new 
observational data, and at the same time develop the basic physics and tools 
for the modeling. The fleld has still a promising future in the third millenium! 
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